arXiv: 1504.06885vl [astro-ph.CO] 26 Apr 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (2015) Printed 28 April 2015 (MN IATeX style file v2.2) 


Searching for Modified Gravity: Scale and Redshift 
Dependent Constraints from Galaxy Peculiar Velocities 

Andrew Johnson^’^’*, Chris Blake\ Jason Dossett^, Jun Koda^, 

David Parkinson^, Shahab Joudaki^ 

^Centre for Astrophysics & Supercomputing, Swinburne University of Technology, P.O. Box 218, Hawthorn, VIC 3122, Australia. 
^ARC Centre of Excellence for All-sky Astrophysics (CAASTRO) 

^INAF - Osservatorio Astronomico di Brera, via Emilio Bianchi 46, 1-23807 Merate, Italy 
^School of Mathematics and Physics, University of Queensland, Brisbane, QLD 4072, Australia 


28 April 2015 


ABSTRACT 

We present measurements of both scale- and time-dependent deviations from the 
standard gravitational field equations. These late-time modifications are introduced 
separately for relativistic and non-relativistic particles, by way of the parameters 
Gma.tter{k, z) and Giight(fc, using two bins in both scale and time, with transition 
wavenumber 0.01 Mpc“^and redshift 1. We emphasize the use of two dynamical probes 
to constrain this set of parameters, galaxy power spectrum multipoles and the direct 
peculiar velocity power spectrum, which probe fluctuations on different scales. The 
multipole measurements are derived from the WiggleZ and BOSS Data Release 11 
CMASS galaxy redshift surveys and the velocity power spectrum is measured from 
the velocity sub-sample of the 6-degree Field Galaxy Survey. We combine with ad¬ 
ditional cosmological probes including baryon acoustic oscillations. Type la SNe, the 
cosmic microwave background (CMB), lensing of the CMB, and the temperature- 
galaxy cross-correlation. Using a Markov Chain Monte Carlo likelihood analysis, we 
find the inferred best-fit parameter values of Gniatter(fc, and Gughtik, z) to be con¬ 
sistent with the standard model at the 95% confidence level. Furthermore, accounting 
for the Alcock-Paczynski effect, we perform joint fits for the expansion history and 
growth index gamma; we measure 7 = 0.665 ±0.0669 (68% C.L) for a fixed expansion 
history, and 7 = 0.73^°;°® (68% C.L) when the expansion history is allowed to deviate 
from ACDM. With a fixed expansion history the inferred value is consistent with GR 
at the 95% C.L; alternatively, a 2a tension is observed when the expansion history is 
not fixed, this tension is worsened by the combination of growth and SNe data. 

Key words: surveys, cosmology: observation, dark energy, cosmological parameters, 
large scale structure of the Universe 


1 INTRODUCTION 

The observation of an accelerating cosmic expansion rate has 
likely provided an essential clue for advancing our theories 
of gravitation and particle physics (Witten 2001). Interpret¬ 
ing and understanding this feature of our Universe will re¬ 
quire both observational and theoretical advancement. Ob- 
servationally it is critical that we both scrutinise the stan¬ 
dard vacuum energy interpretation and thoroughly search 
for unexpected features resulting from exotic physics. Such 
features may exist hidden within the clustering patterns of 
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galaxies, the coherent distortion of distant light rays, and 
the local motion of galaxies; searching for these features is 
the goal we pursue herein. 

Either outcome will facilitate progress: failure to de¬ 
tect unexpected features, confirming a truly constant vac¬ 
uum energy, will give credence to anthropic arguments for¬ 
mulated within String Theory (Susskind 2003). New obser¬ 
vational signatures should then be targeted (e.g., Bousso, 
Harlow & Senatore 2013). Alternatively, an observed devi¬ 
ation from a cosmological constant would indicate a new 
dynamical dark energy component or a modification to Ein¬ 
stein’s field equations (Clifton et al. 2012; Copeland, Sami 
& Tsujikawa 2006). Independent of observational progress. 
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historical trends in science may offer an independent tool 
to predict the fruitfulness of each interpretation (Lahav & 
Massimi 2014). 

The possibility of new physics explaining the accelerat¬ 
ing expansion has inspired an impressive range of alterna¬ 
tive models. As such, a detected deviation from the stan¬ 
dard model will not present a clear direction forwards, that 
is, interpreting such a deviation will be problematic. One 
potential solution, which we adopt, is to analyse observa¬ 
tions within a phenomenological model that captures the 
dynamics of a large range of physical models (e.g.. Bean 
& Tangmatitham 2010; Daniel et al. 2010; Simpson et al. 
2013). It should be noted that not all approaches that in¬ 
troduce modified gravity or dark energy invoke an artificial 
separation between the cosmological constant problem and 
the problem of an accelerating expansion (e.g., Copeland, 
Padilla & Saffin 2012). 

To characterise the usefulness of phenomenological 
models we consider their ability to describe known physi¬ 
cal models: namely, their commensurability (Kuhn 1970). 
This property can be understood as describing the degree 
to which measurements made in one model can be applied 
to others. The absence of this property implies that a mea¬ 
surement should only be interpreted in terms of the adopted 
model: a consistency test. Whereas given this property one 
can constrain a range of models simultaneously, alleviating 
the problem of having to re-analyse each model separately. 

Specifically, the model we adopt allows extensions to the 
standard ACDM model by introducing general time- and 
scale-dependent modifications (Gught and Gmatter) to Gen¬ 
eral Relativity (Daniel et al. 2010): these parameters vary 
the relationship between the metric and density perturba¬ 
tions (i.e, they act as effective gravitational coupling). In 
this case, the equivalence between the spatial and temporal 
metric perturbations is not imposed. The commensurability 
of our model to others can then be shown by proving that 
Giight and Gmatter Capture all the new physics in specihc 
modified gravity scenarios. 

For example, de Felice, Kase & Tsujikawa (2011) show 
that by introducing parameters equivalent to Gught and 
Gmatter oue Can provide an effective description of the entire 
Horndeski class of models. Importantly, the Horndeski class 
of models contains the majority of the viable Dark Energy 
(DE) and modified gravity (MG) models (Silvestri, Pogosian 
& Buniy 2013; Deffayet et al. 2011). An often disregarded 
caveat is that the mappings between these gravitational pa¬ 
rameters and MG and DE theories are only derived at linear 
order. Therefore, until proved otherwise, the ability of the 
phenomenological models to describe physical models is lost 
when using observations influenced by non-linear physics. 
To avoid this reduction in applicability we will focus on 
observations in the linear regime. We note this point has 
been emphasized elsewhere by, for example, Linder & Gahn 
(2007) and Samushia et al. (2014). 

In pursuit of deviations from the standard model we 
use a range of cosmological observations. In particular, two 
dynamical probes will be emphasised: the galaxy multipole 
power spectrum and velocity power spectrum (for example, 
Beutler et al. 2014; Johnson et al. 2014). Hitherto, in the 
context of phenomenological models with scale-dependence, 
neither probe has been analysed self-consistently. In addi¬ 
tion we utilize the following cosmological probes: baryon 


acoustic oscillations, Type la SNe, the cosmic microwave 
background (GMB), lensing of the GMB, and temperature- 
galaxy cross-correlation (this correlation is caused by the 
Integrated Sachs-Wolf effect). 

We adopt this combination of probes, direct peculiar 
velocities (PVs) and redshift-space distortions (RSDs), to 
maximise our sensitivity to a range of length scales. This 
range is extended as the sensitivity of both measurements is 
relatively localised at different length scales: redshift-space 
distortions at small scales, and peculiar velocity measure¬ 
ment at large scales (Dodelson 2003). The benefit is an 
increased sensitivity to scale-dependent modifications. The 
properties of, and physical motivations for, scale-dependent 
modifications to GR are discussed by Silvestri, Pogosian & 
Buniy (2013), and Baker et al. (2014). 

In Section 2 we summarise the adopted phenomenologi¬ 
cal models and further motivate their use. Then in Section 3 
we outline the primary datasets used along with the method¬ 
ology we use to analyse them. Section 4 then presents the 
secondary datasets we employ. The results and interpreta¬ 
tions of the MGMG analysis are presented in Section 5, and 
the conclusions are outlined in Section 6. 

2 MODIFIED GROWTH & EVOLUTION 
2.1 Introduction 

Working within the conformal Newtonian gauge, perturba¬ 
tions to the Robertson-Walker metric can be characterised 
by two scalar potentials. One scalar potential describes a 
temporal perturbation to the metric ip, the other a spatial 
perturbation, (p. The line element in this case is given by 

ds^ = a^[—(1 -I- 2ij))dT^ + (1 “ 2(j))d5j] , (1) 

where a is the scale factor, r is the conformal time - related 
to the proper time of co-moving observers by r = J dt/a{t) - 
and X the spatial coordinate. A non-relativistic fluid within 
this space-time is characterised in terms of a velocity diver¬ 
gence 9(x,t) and a density perturbation Sp{x,T). The cos¬ 
mic evolution of this fluid is then determined by its coupling 
to the metric potentials. 

We concentrate on modifying two of the four gravita¬ 
tional field equations, by requiring energy-momentum con¬ 
servation = 0), or equivalently, by requiring the 

contracted Bianchi identity to hold, i.e., V“G^^ = 0. En¬ 
forcing either constraint one finds the relativistic continuity 
and Euler equations in Fourier space: 

Sm = —9m J- 30, (2) 

dm = —'H9m + k^ip ■ (3) 

where 5m = Spmjpm and H = a/a = (dajdx)!a, and pm is 
the background matter density. This system of four variables 
can then be closed by specifying the gravitational field equa¬ 
tions; in particular, by defining the relationship between the 
two metric potentials, and the coupling between the metric 
potentials and the matter over-density. In GR these rela¬ 
tionships are given by 

= 47rG]va^PmAm (4) 

0 = V', (5) 
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where Gn is Newton’s gravitational constant, and the equa¬ 
tions are defined in terms of the comoving-gauge density 
perturbation = Sm + {3Hlk'^)e m • 

2.2 Glight(^, 2 :) and Gmatter(^, 

We now introduce two dimensionless free parameters Gught 
and Gniatter that we use to model deviations to the field 
equations. Our model is now specified as (Daniel & Linder 
2013) 

= 47rGiva^pmAm X Gmatter ( 6 ) 

V^{(j) + 1p) = SttGnO^ pmAm X Glight ■ (7) 

The first equation governs the motion of non-relativistic par¬ 
ticles, while the second controls the propagation of light 
along null geodesics. As a result, Gmatter can be measured 
using RSDs and direct PVs, and Giight can be measured 
using weak leasing. Because of this distinction the two pa¬ 
rameters are significantly less correlated than models in¬ 
volving a ‘slip’ relation (e.g.. Bean & Tangmatitham 2010). 
Note that the variables {S,/r} in Simpson et al. (2013) and 
Zhao et al. ( 2012 ) are equivalent to {Gught, Gmatter}. There 
is also a trivial re-mapping to the {Q, R} parameters used 
by Bean & Tangmatitham (2010), through Gmatter = QR, 
Giight = <3(1 + 7?)/2. 

To ensure our model can test for a variety of deviations 
from GR we allow for both scale- and redshift-dependence: 
that is, Giight — Giight (^, ^) and Gmatter — Gmatter (^, ^). 
To specify these parameters we use a high vs. low-redshift, 
large vs. small scale binning approach introduced by Daniel 
& Linder (2010). Note, however, that very general func¬ 
tional forms for these parameters (including scale-dependent 
terms) have been developed (Silvestri, Pogosian & Buniy 
2013; Baker et al. 2014). We leave such investigations to 
future work. 

Our adopted model introduces 8 free parameters and 
requires one to specify a redshift and wavenumber transi¬ 
tion scale, zt and kt- We set Zt = 1 and fee = 0.01 Mpc“^; 
therefore, we have two redshift bins (viz., 0 < 2 < 1 and 
1 < z < 2) and two wavenumber bins (10“^Mpc“^ < fe < 
10 “^Mpc“^ and 0.01 Mpc“^ < fe < 0.1Mpc“^), while for 
z > 2 and fe < 10“’^Mpc“^ GR is restored. The transition 
between bins is implemented using an arctan function of 
width Az = 0.05 and Afe = 0.001. 

For our first model we choose to leave the cosmic ex¬ 
pansion unmodified at the ACDM prediction, and concen¬ 
trate on the growth of structure. Henceforth, we will refer 
to this model as model 1. To calculate the relevant observ¬ 
ables (to be discussed in the next section) we use camb and 
CosmoMC. The modified field equations (Eq 7) are incorpo¬ 
rated into camb using the publicly available code ISITGR 
(Dossett, Ishak & Moldenhauer 2011), and the exact equa¬ 
tions implemented in camb are given by Dossett, Ishak & 
Moldenhauer (2011). Note the only significant difference be¬ 
tween the equations employed in camb and Eq (7) is that 
the latter are written within the synchronous gauge (Ma & 
Bertschinger 1995). 

A few technical comments on the model are unavoid¬ 
able: Firstly, super-horizon curvature perturbations need to 
be conserved independent of the form of field equations 
(Bertschinger & Zukin 2008). This condition was shown to 


be satisfied for this model by Pogosian et al. (2010). Addi¬ 
tionally, it is natural to include a smoothness theory prior 
on these parameters, however, given the large distance be¬ 
tween the centre of our bins we choose not to include such a 
prior (Silvestri, Pogosian & Buniy 2013). With more accu¬ 
rate data, and hence a larger number of bins, this argument 
will no longer be valid. Finally, the accuracy of any map¬ 
ping from our model to physical models (i.e., those derived 
from an action) relies on the validity of the quasi-static ap¬ 
proximation (QSA). Following the arguments presented in 
Silvestri, Pogosian & Buniy (2013) it is reasonable to include 
a theoretical prior to ignore such deviations. 


2.3 Varying Growth and Expansion: {7,wo,uia} 

As more freedom is introduced to model deviations from 
GR the precision of the inferred parameters degrades. We 
must decide then which features of the standard model to 
preserve; for example, to what extent does the expansion 
history dictate the growth history. This presents a balancing 
problem with no clear solution. To partially circumvent this 
issue we adopt a second model (which we label model II) . In 
contrast to our first model, this model includes only minimal 
extensions to the standard model. As a result there are fewer 
free parameters and more precise tests are possible (although 
we nonetheless introduce deviations to both the expansion 
and growth history). 

This minimal extension to the standard model using 
the parameters {uio,uia, 7 } has been advocated by Linder 
& Gahn (2007); Linder (2005), and Simpson & Peacock 
( 2010 ), and applications have been presented, for exam¬ 
ple, by Huterer & Linder (2007). To expand on this, we 
introduce deviations to the expansion history through a 
time-dependent equation of state w{z), which is expressed 
in terms of two free parameters: wo = u>(a = 0 ) and 
Wa = —(dw/da)l^_^, as a function of the redshift w(z) = 
wo+Waz/(l + z). Note the expansion history is still governed 
by the Friedman equation, there is simply more freedom 
in the properties of the dark energy component. We intro¬ 
duce deviations in the growth history by parameterizing the 
growth rate as f{z) = i2rn{z)'^, where 7 is the growth index; 
within GR one expects 7 ~ 0.55. The growth rate is defined 
by /(a) = d In D(a)/ci In a, and D{a) = 8{a)/8{a = 1). 


3 PRIMARY DATASETS: METHODOLOGY 

Below we will outline the measurements we use in Sec. 5, 
in addition to the tools we use to analyze them. A gen¬ 
eral summary is provided in Table 1 where the datasets, 
the measured quantities, and the fitting ranges adopted are 
specified. The focus will be on introducing extensions to the 
public MGMG code CosmoMC (Lewis & Bridle 2002) and camb 
(Lewis, Ghallinor & Lasenby 2000) to update the range of 
datasets one can analyze. 


3.1 Velocity Power Spectrum 

The radial PVs of galaxies in the local universe induce a 
fluctuation in the apparent magnitude m, defined as (Hui & 
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Table 1. Summary of the datasets used in this analysis. Given model I includes scale-dependent terms, we divide our measurements into 
three separate groups: those used to constrain model I &; II, only model I, and only model II. This division is indicated by the horizontal 
lines, and follows the order in which the categories were introduced. 


Cosmological Probe 

Dataset 

Measured quantity 

Reference 

CMB temperature. . . 

Planck 

CfT 

Planck Collaboration et al. (2013a) 

CMB polarization... 

WMAP-9 

E 

Bennett et al. (2013) 

CMB-Lensing. 

Planck 

(j4> 0 

Planck Collaboration et al. (2013b) 

BAOs. 

6dFGS 

rsIDviz) 

Beutler et al. (2011) 


BOSS DRll LOWZ 

Dv(r“/rs) 

Anderson et al. (2013) 


BOSS DRll QSA-Lya 

H(z)rs, Da/ts 

Font-Ribera et al. (2014) 


BOSS DRll Lya 

H(z)rs,DA/rs 

Delubac et al. (2014) 

Type la Supernovae .. . 

SNLS 

fi{z} 

Conley et al. (2011) 

Dataset extension I 

ISW-density cross. .. 

WMAP3 


Ho et al. (2008) 

Velocity Power Spectrum 

6dFGSv 

’Pvvik) 

Johnson et al. (2014) 

BAO (reconstructed)** 

WiggleZ 

Pv(rflyr.) 

Kazin et al. (2014) 


DRll CMASS 

DA(^)(r“/r,),H(z)(r,/rfld) 

Anderson et al. (2014) 

Power Spectrum Multipoles DRll CMASS 

Po{k),P2{k) 

Beutler et al. (2014) 


WiggleZ {zef[ = 0.44) 

Po{k),P2ik),P4{k) 

Blake et al. (2011a)^ 


WiggleZ (zeff = 0.73) 

Po(fc),P2(fc),P4(fc) 

Blake et al. (2011a)^ 

Dataset extension II 

RSDs . 

6dFGS 

/o-8(z) 

Beutler et al. (2012) 

RSD-BAO-AP . . . 

WiggleZ 

A(z),FAp{z),fas{z) 

Blake et al. (2012) 

RSD-BAO-AP . . . 

BOSS CMASS 

■Dv As ( 2 ), Pap( 2 ), frrsiz) 

Beutler et al. (2013) 


“ Both the reconstructed BAO measurements (CMASS and WiggleZ) have been calculated by marginalising over the general shape of 
the correlation function. Marginalising over the shape decorrelates the BAO measurement with the power spectrum multipole 
measurement, allowing one to fit for both measurements simultaneously. 

Note, however, these measurement have been updated in this work using an improved methodology. 


Greene 2006) 

5m(z) = lm(z) — rh{z)] . ( 8 ) 

The over-bar indicates that the variable is being evalnated 
within a homogeneous nniverse, namely, a universe with no 
density gradients and therefore no pecnliar velocities. Recall 
the apparent magnitude is defined as 

m = M-b 51ogiQ(DL(2)) + 25 . (9) 

Here M is the absolute magnitude, and Dl(«) the luminos¬ 
ity distance. The presence of large scale clustering induces 
fluctuations in Sm{z) from galaxy to galaxy (this is equiva¬ 
lent to a peculiar velocity), furthermore, these fluctuations 
are correlated for nearby galaxies (Hui & Greene 2006; Gor¬ 
don, Land & Slosar 2007). The magnitude of both effects 
can be described by a covariance matrix which we define 
as Cfj = {Smi(zi)Smj(zj)} . Once a model is specihed this 
covariance matrix can be calculated as 

C^ = Gizi,Zj) I ^k^V..{k,a^l)W{k,aij,n,ri). 

( 10 ) 

Where Vvvik) — 'Pee{k)/k^ is the velocity power spectrum, 
and 6 = V • u is the velocity divergence, furthermore 

W{k,aij,ri,rj) = l/3[jo{kAij) - 2j2ikAij)] fi-fj 

1 2 

+ Aij)rirj sin {atj) , 

G{zi,Zj) = 

( 5 y / _ {i+zj^ \ ( _ {i+z,f \ 

Vlnio; V H{zJDUz^)J V H{zj)DUzj)J ’ 


where aij = cos~^{fi ■ fj), Aij = |ri — r^j and is the po¬ 
sition vector of the i**' galaxy. This analytic solution for the 
window function was presented by Ma, Gordon & Feldman 
(2011). For further details on this calculation we refer the 
reader to Johnson et al. (2014). 

We perform a full likelihood calculation using the 
6 dFGSv peculiar velocity sample (Springob et al. 2014). To 
calculate the covariance matrix in Eq. (10) we integrate over 
the wavenumber range k = 0.0005 — 0.15/i Mpc“^. Given 
the dominance of large-scale information in peculiar velocity 
measurements, we neglect velocity bias in this calculation. 

In order to minimise the influence of poorly understood 
non-linear effects a non-linear velocity dispersion component 
crpv is introduced into the diagonal elements of the covari¬ 
ance matrix (Silberman et al. 2001). This nuisance parame¬ 
ter is marginalised over in the analysis. The covariance ma¬ 
trix is thereby updated: 

Eij = Gij -b tTpySij . (11) 

One can now define the posterior distribution as z 

P(E|5m) = |27rE|~^^^ exp , (12) 

where Jm is a vector of the observed apparent magni¬ 
tude fluctuations. Note the dependence on the cosmological 
model is introduced through the covariance matrix. 

The model velocity power spectrum is generated using 
a transfer function. This can be defined starting from the 
peculiar velocity in the synchronous gauge Vp^ (cf., Ma & 
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Bertschinger 1995)^. As this gauge is defined in the dark 
matter rest frame, i.e., there are no temporal goo perturba¬ 
tions, a gauge transformation is necessary. Using the con¬ 
vention of Ma & Bertschinger (1995) we define h and rj as 
the metric perturbation in the synchronous gauge. Now by 
moving into the Newtonian gauge one finds the appropriate 
transfer function: 

Tv{k) = -G(^ka +pbvfi'>l{pb+, (13) 

where k^a = h/2 -\- Sp. 

In Fig. 1 we plot the measurements of Vvvik) by John¬ 
son et al. (2014), here the blue (green) points were mea¬ 
sured using the OdFGSv (low-2; SNe) sample. For this plot 
the black line shows the power spectrum prediction assum¬ 
ing GR, while the red and orange lines show the predictions 
for different values of the post-GR parameters. For these cal¬ 
culations the Planck best-fit parameters are assumed. Ad¬ 
ditionally, the green line shows the prediction when using 
our best-fit parameter values (see sect. 5 for details). Note 
the time evolution of the density perturbation Am is set 
by a friction term 2'HAm and a source term There¬ 

fore, by modifying Gmatter one changes the source term to 

k'^'lp ~ a^Gmatter(fe, 2;)Am; heUCe, with Gmatter (fc, 2:) > 1 

both the late-time clustering and the amplitude of the ve¬ 
locity power spectrum are enhanced. 

3.2 Power Spectrum Multipoles 

We measured the multipole power spectra of the WiggleZ 
Survey data using the direct estimation method introduced 
by Yamamoto et al. (2006) and extended by Blake et al. 
(2011a) and Beutler et al. (2014). We provide a brief sum¬ 
mary of the technique here, referring the reader to the above 
papers for a full description. 

The redshift-space 2D galaxy power spectrum Pg{k,p), 
where p is the cosine of the angle of the wavevector k with 
respect to the line-of-sight, may be expressed in terms of 
multipole moments Pe{k) using a basis of Legendre polyno¬ 
mials Le{p): 

p;(fc,/i)= Pr(fc)Lr(/i), (14) 

even i 

where 

2 / 4-1 

Pi{k) = —-— J dpPfik,p) Le{p). (15) 

The power spectrum multipoles provide a form of data com¬ 
pression; in linear theory all the information is contained in 
the £ = 0,2,4 terms, with the first two multipoles dominat¬ 
ing the observed signal. 

The rapid estimation technique of using Fast Fourier 
Transform (FFT) methods to measure Pfik,p) in bins of k 
and /i, where p is defined with respect to a fixed axis parallel 
to the line-of-sight of the field centre, and then estimating 
Pi{k) by a direct sum over the binned results using Equa¬ 
tion 15, has two difficulties. First, for a wide-area survey 
the line-of-sight direction with respect to which p should be 
measured will not be fixed. Secondly, at low k the sum over 
p bins is problematic to evaluate due to the limited number 

^ Our starting point is set by variables used within camb. 


of modes available in Fourier space. The Yamamoto et al. 
(2006) method estimates Pi{k) using a sum over all galaxies 
for each wavevector k on the FFT grid, allowing the line- 
of-sight vector to vary for each object and without binning 
in p. Window function effects are included using a similar 
sum over unclustered objects. Additive corrections are in¬ 
cluded for shot noise and for the discreteness of the grid. 
The measurements are then binned by wavenumber k = |fc|. 

Following the analysis of the WiggleZ baryon acoustic 
oscillations (Blake et al. 2011b), we estimated the £ = 0, 2,4 
multipole power spectra in the (9,11,15, 22,1, 3)-hr survey 
regions in the overlapping redshift ranges 0.2 < 2: < 0.6, 
0.4 < 2: < 0.8 and 0.6 < 2 < 1.0. We measured the spec¬ 
tra in 14 wavenumber bins of width Ak = 0.02 ft Mpc“^ in 
the range 0.02 < k < 0.3ft Mpc“^. For this analysis, how¬ 
ever, we only use the non-overlapping redshift ranges that 
we label low-z and high-z. The results for the monopole and 
quadrupole are given in Fig. 2. 

We determined the covariance matrix of each vector 
\Po{k), P 2 {k), Piik)] by repeating the measurements in each 
survey region for a series of 600 mock catalogues, built from 
N-body simulations generated by the method of COmov- 
ing Lagrangian Acceleration (GOLA; Tassev, Zaldarriaga & 
Eisenstein 2013). As described by Kazin et al. (2014) we 
produced a halo catalogue by applying a friends-of-friends 
algorithm to the dark matter particles, and populated the 
haloes with mock galaxies using a Halo Occupation Distribu¬ 
tion such that the projected clustering matched that of the 
WiggleZ galaxies. The mocks were sub-sampled using the 
selection function of each region, and galaxy co-ordinates 
converted to redshift-space. 

We also determined the convolution matrix for each 
region and redshift slice, which should be used to project 
a model multipole vector to form a comparison with the 
data given the survey window function. For a wide-angle 
survey such as the BOSS, determination of the convolution 
involves a numerically-intensive double sum over randomly- 
distributed objects (Beutler et al. 2014). However, for the 
more compact WiggleZ Survey geometry, we found that it 
was acceptable (in the sense that any offset was far smaller 
than the statistical error) to use a flat-sky approximation, in 
which FFT methods were used to convolve a series of unit 
multipole vectors, generating each row of the convolution 
matrix in turn. 

In addition to the WiggleZ multipole measurements, we 
include the monopole and quadrupole measurements from 
the BOSS-DRll CM ASS sample presented in Beutler et al. 
(2014); the reader is referred to this paper for technical 
details on the calculation. From the CMASS sample the 
1 = 0, 2 multipole power spectrum are calculated for the 
wavenumber range k — 0.01 — 0.20ft Mpc“^with a spacing 
of Ak = 5 X 10“^ft Mpc“^. These measurement are pre¬ 
sented for both the North and South Galactic Cap regions 
at an effective redshift of 2eff = 0.57. 

We plot the CMASS multipole measurements in Fig. 3. 
For this plot the blue-dashed (red-dashed) lines show the 
multipole predictions when setting Gmatter (ft > 0.01; 2 < 
1) = 1.8 (G matter (ft > 0.01; 2 < 1) = 0.3), while the black 
lines show the prediction assuming GR. For these predictions 
the best-fit parameters from Planck are assumed, in addition 
we set the bias to 6 = 1.85, shot noise to N = 1800 ft“®Mpc®, 
and the velocity dispersion to cr„ = 4ft“^Mpc. Moreover, 
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Figure 1. The velocity power spectrum Vwik) at 2 = 0 for different parameter combinations of the adopted phenomenological model. 
The black line shows the prediction assuming General Relativity, and the orange and red lines illustrate the effect of varying the low -2 
and high-/c bin for Gmatter- For the red line Gmatter(' 2 : < 1; fc > 0.01) = 1.8 and for the orange line Gmatter(^ < 1; /c > 0.01) = 0.3: for 
these predictions the standard cosmological parameters are fixed at the Planck best-fit values, and unless specified otherwise all non-GR 
parameters are set to be consistent with GR (i.e., set equal to 1). Moreover, the green line shows the prediction found using the best-fit 
parameter values found using set 4 (see sect. 5 for details). The best-fit values here correspond to the parameter values that maximise the 
likelihood. The blue and green data points correspond to the 68% confidence intervals for the mean power within each bin for the GdFGSv 
data and the low -2 SNe data set constructed in Johnson et al. (2014). The thick black line indicates the the mean power predicted by 
GR in each k-bin, this is calculated assuming a Planck cosmology. 


the orange lines give the prediction when using our best-fit 
model parameters (see sect. 5 for details) Note, for simplicity 
the theory predictions have only been convolved with the 
NGC window function. 


3.2.1 Modelling the Power Spectrum Multipoles 

To model the redshift-space 2D galaxy power spectrum 
Pg{k,fi) we use linear theory plus an empirical Gaussian 
damping term (Hatton &: Cole 1998); the resulting model is 
given by 

Pg{k,f,) = [Pgg{k) - 2^ VgS (fc) + /^69 (fc)] D{n,k) , 

(16) 

where D{ij,,k) = exp[—(fc/^Cv)^]. The standard interpreta¬ 
tion of this damping, which is clearly observed in redshift 
surveys, is the uncorrelated pairwise velocity dispersion of 
galaxies. We absorb our ignorance by treating CTv as a free 
parameter to be marginalised over for each survey. 

Assuming linear theory the continuity equation (eq. 3) 
can be written in Fourier space as 


e{k) = -f{a)S{k). 


(17) 


However, we are modifying the gravitational held equa¬ 
tions, so one needs to be self-consistent, given that the 
modihcations (Eq. 7) will change the growth rate in a 
scale-dependent manner. We calculate this modihed scale- 
dependent growth rate as 


f{k,a) 


din Ac(fe, a) 
din a 


(18) 


This is self-consistent given camb contains all the relevant 


physics, i.e, the density and velocity variables are evolved 
according to the modihed held equations. As a reminder of 
the potential scale-dependence we write the growth rate as 
f{k). Now assuming a local, scale-independent linear bias 
(dg = bS) and no velocity bias {Og = 9) Eq. (16) reduces to 

P^{k, m) = (Pssik) + N){1 + f{k)^^/hf 7?(/i, k) , (19) 

Here we have included a shot noise component N, this 
is treated as a free parameter for the CMASS analysis. 
To justify the previous assumptions we truncate the ht 
for both the WiggleZ and CMASS multipoles to relatively 
large scales; to wit, we set = O.lOh Mpc“^and 

= O.lS/i Mpc“^. The WiggleZ measurements are 
used to a higher wavenumber because of the smaller bias of 
the sample (6 ~ 1), in addition to the larger error bars^. 
The matter power spectrum is calculated within camb using 
only linear theory: we choose not to incorporate non-linear 
corrections via HALOFIT. The use of HALOFIT presents an is¬ 
sue as the corrections have not been shown to be valid for 
general modihed gravity models. 

In order to correctly interpret RSD measurements one 
is required to consistently incorporate our ignorance of the 
expansion history of the universe (viz., H{z)), bearing in 
mind that these measurements are performed assuming a 
hducial cosmological model. As a result, in a trial cosmol¬ 
ogy, the growth rate measurements should be adapted us¬ 
ing the covariance with the Alcock-Paczynski (AP) distor¬ 
tion. Any discrepancy between the chosen hducial expan- 


^ With a lower biased tracer, for example, the effect of non-local 
halo bias is less significant (Chan, Scoccimarro &; Sheth 2012). 
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Figure 2. The monopole, quadrupole and hexadecapole measurements from the WiggleZ survey, for both the high- 2 : and low -2 samples 
(2 = 0.44,0.73, respectively). For simplicity we combine the results from the 6 different survey regions; however, note this is not the 
format of the data we use: each survey region has a different window function and therefore is analysed separately. 


sion history {Da(z), H{z)) and the physical expansion his¬ 
tory (Da(z), H{z)) can be accounted for by scaling the true 
(physical) radial and tangential wavenumbers 
The amplitude of the wavenumber scalings is determined by 


ay 


Da{z) 

H{z) ’ Df{z) 


( 20 ) 


Hence the observed wavenumbers are given by = 

and Including this scaling in Eq 

(19) one finds (Ballinger, Peacock & Heavens 1996; Matsub- 


ara & Suto 1996; Simpson & Peacock 2010) 


P|(fc',/r') = 


1 I '■ 

1 + /r 


1 fj, 


/2 / 01 


1 + /? 
a\la\ 
1 -2 


- 1 


- 1 


X Pss 


k' 


1 + n'- 


- 1 


( 21 ) 

X -D(/r, k) 


where k' = + (^J|*’°)^i fi ~ a^nd P = f/b. 

This scaling introduces a new source of anisotropy in the 
clustering of galaxies, making it partially degenerate with 
redshift-space distortion effects, accordingly it is important 
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to account for this effect in this type of analysis (Blake et al. 
2012; Beutler et al. 2014). 

Two components must be included to compare our the¬ 
oretical predictions with observations: the window function 
and integral constraint effect, both of which result in a dis¬ 
tortion to the measured power spectrum relative to the true 
power spectrum. Window function effects are induced by the 
complex geometry of the survey (viz, a non-cubical geome¬ 
try); and the integral constraint effect occurs as the condi¬ 
tion 5k=o = 0 is applied to the data: this imposed normal¬ 
ization for the k — 0 mode is invalidated by super-survey 
modes. Both effects induce a suppression of power at low-fc 
(Peacock & Nicholson 1991; Beutler et al. 2014). 

A consistent comparison between our model and the 
observations therefore requires us to include the window 
function effects in our modelling. Following Beutler et al. 
(2014) the convolved multipoles Pf°'"'(fc) are calculated for 
the CMASS sample as 

prjk) = 27v f 

L 

( 22 ) 

where 

iVran 

\W{k, k )\ii^ = 2i (—i) (2£-f 1) ^ ( WFKp(a:i)®FKp(a:j) 

jt{k\^:^x\)jL{k'\^x\)Ct{xh ■ Ax)CLixh ■ Ax), 
and the integral constraint term is given by 


this estimator biased (Hartlap, Simon & Schneider 2007). 
We correct this bias using the estimator 


c-^ 


Ns - rib-2 
Ns - 1 


r* 

'^mock 


(23) 


where rib is the number of power spectrum bins, and Ns the 
number of mock realisations used to construct the covariance 
matrix. 


3.3 BAOs 

Acoustic oscillations in the photon-baryon plasma, prior to 
recombination, imprint a series of fluctuations in large-scale 
structure: in configuration-space one finds a preference for 
galaxies to be distributed with a given comoving separation 
(~ 105 /i“^Mpc). 

This excess in clustering (the BAO feature) functions as 
a cosmic yard-stick allowing the cosmic expansion history to 
be mapped out. By measuring the spherically averaged BAO 
position one determines 

Dv{z) = [cz{l + z)^Da{z)^H{z)] . (24) 

Here Da{z) is the angular diameter distance. With higher 
signal to noise measurements one can extract more infor¬ 
mation by isolating the transverse and line-of-sight BAO 
positions, determining 

ctperp = DA{z)rs‘^/DA^{z)rs (25) 

apar = Pfld(z)rf/P(2)rs. (26) 


Pf(fc) = 2^ 


\Wik)\j 

|W(0)|g 




Here Jl are spherical Bessel functions of order L, Wan is the 
number of galaxies in the synthetic catalogue, and we sum 
over the monopole and quadrupole {L = 0, 2). 

Each survey region has a different window function 
and hence needs to be treated separately. To compute the 
CMASS likelihood we use the publicly available CMASS 
window functions^. The combined likelihood is now com¬ 
puted as 


-21n(r'^'‘5giez^ - 21n(£®°®®) = 

^(^Wigglez _ 
i=l 

2 

E(^ 


3BOSS _ pConv'jT^-1 


BOSS,j\^j 


- p. 


^pBOSS pConv^ 


The i indices specify the two redshift bins {zes = 
0.44,0.73) and six survey regions for WiggleZ (12 sepa¬ 
rate measurements). The j indices specify the two sur¬ 
vey regions (NGC and SGC) for CMASS. Furthermore, 

^WiggleZ ^ pconv^^^]. pSOSS ^ 

[Po"°"''(fc),P|°"'-(fc]]j. 

The hat in C ^ indicates that we are using a statistical 
estimator for the inverse covariance matrix. This estimator 
is determined by the covariance matrix measured from mock 
catalogues: typically one would use C~^ = however, 

the noise in the derived covariance matrix (C'”^^,]^) makes 


^ https://sdss 3 .org/science/boss_publlcations.php 


By including the dependence on Cg (the sound horizon 
at the drag epoch), and expressing the measured quantity 
as a ratio of the fiducial prediction, the dependence on CMB 
physics and the assumed cosmology has been made explicit. 

To constrain the expansion history we use the follow¬ 
ing BAO measurements: WiggleZ reconstructed from Kazin 
et al. (2014), reconstructed DRll-CMASS and DRll- 
LOWZ from Anderson et al. (2014), and the 6dFGS mea¬ 
surement from Beutler et al. (2011). By ‘reconstructed’ we 
are referring to the process of sharpening the acoustic peak 
by using information from the local density field (cf. Pad- 
manabhan et al. 2012). The above measurements (excluding 
CMASS) can be incorporated into a likelihood given by 

-2lnC = {x-S)'^C~4x-S), (27) 

with the theory vector 

X = [Dv(0.44)(rfid/ra), Dv(0.6)(rfid/rs), I>v(0.73)(rfid/?'s) 

T»v(0.32)/rd,rg/T>v(0.106)], (28) 

the data vector 


5 = [1716, 2221, 2516,8.25, 0.336], (29) 


and the covariance 

matrix^ 





/2.17898 

1.11633 

0.46982 

0 



1.11633 

1.70712 

0.71847 

0 

0 

^-1 _ 
^BAO — 

0.46982 

0.71847 

1.65283 

0 

0 


0 

0 

0 

36.025 

0 


( 0 

0 

0 

0 

4444.4/ 


We have scaled the WiggleZ elements for clarity; the true co- 
variance matrix is obtained by scaling the WiggleZ elements by 
10-“: = 2.17898 X 10"“. 
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Figure 3. The monopole and quadrupole power spectrum for both the BOSS-DRll CMASS survey regions (NGC and SGC). The 
blue-dashed line shows the prediction with Gmatter(fc > 0.01; 2 : > 1) = 1.8, the red-dashed lines Gmatter(fe > 0.01; 2 < 1) = 0.3; for 
these predictions the best-fit parameters from Planck are assumed and we set the bias to fe = 1.85, shot noise to N = 1800 h“®Mpc®, 
and the velocity dispersion to (t„ = 4h~^Mpc. For simplicity these theory predictions have only been convolved with the NGC window 
function. The orange lines gives the prediction from the best-fit model parameters (see sect. 5 for details), convolved with the NGC 
window function. Note, for the final analysis we only fit our model to fcmax = O.lO/i Mpc“'^. 


The CMASS measurements are in the form of probabil¬ 
ity distributions for P(Q:perp) and P(Qpar) evaluated at 
2eff = 0.57. These measurements are therefore analysed sep¬ 
arately, for details see Anderson et al. (2014). A number of 
these BAO measurement have been calculated using the ap¬ 
proximate fitting formula for rs(2d) from Eisenstein & Hu 
(1998); hence throughout, where appropriate, the BAO mea¬ 
surements derived using this approximation are scaled to be 
consistent with the result from camb (cf., Mehta et al. 2012). 

To further improve the redshift range of our expan¬ 
sion history measurements we extend this ‘base’ sample by 
including the Lyman-a BAO measurement from Delubac 
et al. (2014), and the Quasar-Lya cross-correlation measure¬ 
ment from Font-Ribera et al. (2014). The measurements are 
Dh{z = 2.34)/r, = 9.18 ± 0.28, Da{z = 2M)/n = 11.28 ± 
0.65, PH(a = 2.36)/ra = 9.0 ± 0.3, T)a(« = 2.36)/rs = 
10.8 ±0.4 , where Du = c/H. Both common cosmic variance 
or a common source for the measurement error would in¬ 
duce correlations between the Lyman-a measurement. For¬ 
tunately, the origin of the dominant error components for 
these measurements are distinct, and hence the measure¬ 
ments are uncorrelated (Font-Ribera et al. 2014). Addition¬ 
ally, we treat any correlations between the BOSS and Wig- 


gleZ surveys as insignificant, given the small overlapping 
area (~ 550deg^) and the significance of shot noise in Wig- 
gleZ measurements. 


3.4 Growth Rate and Alcock-Paczynski 
Measurements 

The growth rate measurements presented in this section will 
be used to constrain 7. Following the arguments presented 
in subsection 3.2.1 we only include growth rate constraints 
that have consistently incorporated the Alcock-Paczynski 
effect. The exception to this point is for very low-redshift 
observations, which are effectively insensitive to changes in 
the expansion history. 

In order to self-consistently express the degeneracy with 
the expansion history we chose to fit to joint 3D posterior 
distributions from AP, BAO and RSD measurements: as op¬ 
posed to marginalized ID constraints on fas{z). The growth 
rate measurements we utilize are measured from BOSS- 
DRll survey, the WiggleZ Dark Energy Survey, and the 
6dF Galaxy survey (Beutler et al. 2014; Blake et al. 2012; 
Beutler et al. 2012). For the CMASS sample we use the data 
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vector® 

Sk°lU.2o = pv(0.57)/r4«d),-FAp(0.57),/(0.57)a8(0.57)] 

= [13.88,0.683,0.422], (30) 

Where the AP effect translates into a geometric constraint 
on Fap{z) = (1 + z)Da(z)H{z)/c. And the corresponding 
covariance matrix is given by 

/36.400 -2.0636 -1.8398\ 

10"C7r°r=o.2o = 0 1.0773 1.1755 (31) 

\ 0 0 2.0438 / 

The WiggleZ survey measurements are performed within 
three overlapping, hence correlated, redshift bins at ZeS = 
0.44, 0.60,0.73. We first split the data vector into redshift 
bins, namely = {Szi,Sz 2 ,Sz 3 ). In each of these 

redshift bins Blake et al. (2012) measure the parameter com¬ 
bination 

= [A{zi), FAp{zi), f{zi)as{zi)] , (32) 

where A{z), the acoustic parameter, is given by 
A{z) = 100^v(z)V»m h2 

The measured values are now = (0.474,0.482,0.413), 
5,^2 = (0.442,0.650,0.390), and 5z3 = (0.424,0.865,0.437). 
Table 2 in Blake et al. (2012) gives the full covariance matrix 
for 5WiggieZ 

The final measurement we use is f(0A06)as(0A06) — 
0.423 ±0.55 from Beutler et al. (2012). As noted previously, 
the AP effect is not significant for this measurement given 
the low-redshift nature of the sample. All of the introduced 
measurements are now incorporated using the likelihood 

-21n£ ^ (x-S)'^C~^(x-S), (34) 

here x, S and C are the appropriate theory vector, data 
vector and covariance matrix. Note that BAO information 
is included in both Sec. 3.3 and Sec. 3.4 and we do not 
double-count this information. 


4 SECONDARY DATASETS 

A brief introduction and motivation is given for the addi¬ 
tional datasets we use. 

4.1 Type-la SNe 

Sample variance effectively imposes a minimum volume limit 
for BAO detection. Accordingly, large volumes and hence 
higher redshift observations are preferable. Type-la SNe 
measurements do not have this restriction and hence can 
provide very accurate constraints on the low-redshift expan¬ 
sion rate: an epoch where the presence of “dark energy” 
appears to dominate. 

Therefore we include the distance modulus measure¬ 
ments for 473 type la SNe presented in Conley et al. (2011). 
The ” SNLS” sample is a combination of a number of previ¬ 
ous surveys combining supernova legacy survey results with 

® This result is found fitting the power spectrum multipoles to 
fcmax = 0.20 h Mpc“^. 


other low-2 and high-2 observations. These measurements 
are included in our analysis using the cosmomc likelihood 
module provided by Conley et al. (2011)®. This likelihood is 
evaluated by (firstly) calculating the model apparent magni¬ 
tudes (or more accurately, the rest-frame peak B-band mag¬ 
nitude) : 

m-rnodei = 51ogior>L(«CMB, «Hei, • • •) “ a(5 - 1) + fiC + Mb 

Here TOp is luminosity distance with the dependence on the 
Hubble constant removed (it’s dimensionless). And 2cmb 
and 2Hei are the CMB frame and heliocentric frame redshifts 
of the SN. A4b is a parameter which controls the zero-point 
and is a function of both the absolute magnitude of the SN 
and Hq, this parameter is marginalised over. The brightness 
of each SN is ‘standardised’ using observations of the shape 
of the light curve, s, and the colour C; in addition to the 
empirical relationship of these parameter with the luminos¬ 
ity of the object: these dependences are characterised by the 
parameters a and /3. 

Writing the model predictions as a vector rn-modei the 
likelihood is given by 

2 In T = (rUobs ^model) C (rU-obs rUixiodel) , (35) 

where fhobs is a vector of the observed B-band magnitudes. 
The elements of the non-diagonal covariance matrix C in¬ 
cludes contributions from the following effects: the intrinsic- 
scatter of type la SN, the errors on the fitted light curve 
parameters, the redshift error, a host correction error, and 
the covariance between s, C and mobs- There are additional 
corrections for the local peculiar velocity field, for further 
details see Conley et al. (2011). 

In Section 5.3 we adopt a second SNe dataset, namely 
the JLA sample (Betoule et al. 2014). This sample is com¬ 
posed of recalibrated SN la light-curves and distances for 
the SDSS-H and SNLS samples; this sample can be distin¬ 
guished from the SNLS sample by the treatment of system¬ 
atic affects, the end result is a 1.8cr shift from the SNLS 
3-year results. 

4.2 CMB 

For the models we adopt GR is restored at the time of the 
last scattering surface; accordingly, the components of the 
temperature fluctuations, unmodified by large-scale struc¬ 
ture, provide a powerful tool to both constrain the physical 
components of the universe and the initial conditions which 
seed large-scale structure. 

The likelihood code for the power spectrum from 
Planck is a hybrid: it is divided into high-f and low-L For 
high-/ (I > 50) we use the likelihood code CamSpec described 
by Planck Collaboration et al. (2014). This algorithm uses 
temperature maps derived at 100, 143 and 217 GHz. Once 
both diffuse Galactic emission and Galactic dust emission 
are masked, 57.8% of the sky remains for the 100 GHz 
map and 37.3% for the remaining maps. At low multipoles 
(2 < I < 49) the likelihood is computed using the CommEuider 
algorithm (Eriksen et al. 2008) using the frequency range 
30-353 GHz over 91% of the sky. 

Sub-Hubble modes near reionization are damped by 

® https://tspace.library.utoronto.ca/handle/1807/25390 
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Thomson scattering, thus obscuring our view of the primor¬ 
dial power spectrum: We observe a fluctuation amplitude 
The degeneracy between the optical depth r and 
the amplitude of the primordial power spectrum As can be 
partially broken by including polarization data: the relative 
amplitude of the polarization and temperature power spec¬ 
trum constrain r. For this purpose we include the large-scale 
polarization measurements (Cf"^) from WMAP-9 (Bennett 
et al. 2013). We use the likelihood code from Planck which 
fits to the Z-range (2 < Z < 32). 

4-2.1 CMB Lensing 

Photons travelling from the last scattering surface to our 
satellites encounter a number of over- and under-densities 
along the way. The intersected structure deflects the photon 
paths and the large-scale clustering of matter causes these 
deflection paths to be correlated over the sky (Blanchard & 
Schneider 1987). The combined effect of this CMB lensing 
is a re-mapping of the CMB temperature fluctuations (cf., 
Lewis & Challinor 2006): 

T(n) = -b V<l>(h)) . (36) 

Where "I?)!!) is the CMB lensing potential given by 

'^'(n) = - / dx G(x, X*) Vo - x) + Six^', Vo - x)] ■ 
Jo 

Here x is the conformal distance, rj is the conformal time (rjo 
is the time today), and G(x, X*) is a weighting function. The 
integration is taken from the last scattering surface (x*) to 
today (x = 0)i hence this term represents the integrated ef¬ 
fect of structure on photon paths, or more accurately, since 
we are interested in testing GR, the integrated effect of spa¬ 
tial and curvature perturbations. 

The lensing power spectrum Cf can be extracted from 
CMB maps; here we use the results from Planck Collabora¬ 
tion et al. (2013a) for the Z-range 40 < Z < 400 (with the 
bin size AZ = 64): this Z-range is chosen as it encompasses 
the majority of the lensing signal (~ 90%) and is likely less 
influenced by systematic effects (cf., Planck Collaboration 
et al. 2013a). Given the lensing kernel peaks at 2: ~ 2 and 
we are only using Z < 400, the lensing power spectrum mea¬ 
surements used are only probing linear scales. Accordingly, 
we use linear theory to predict the lensing power spectrum 
and expect no systematic errors to be introduced from this 
modelling. 

4-2.2 Temperature-Galaxy Cross-Correlation 

At late times the accelerating cosmic expansion dictates the 
evolution of density perturbations, one consequence is time- 
dependent metric potentials. This time-dependence is ap¬ 
parent in the CMB as it generates a net energy loss for 
CMB photons as they propagate through these potential 
wells (Sachs & Wolfe 1967). This feature is known as the 
integrated Sachs-Wolfe (ISW) effect. The influence on the 
CMB power spectrum is given by 

+ (38) 

The ISW effect induces a correlation between the CMB 
(low-Z) and large-scale structure probes: this is measured 


using the temperature-galaxy cross-correlation power spec¬ 
trum Gf^ (cf. Ho et al. 2008). For our analysis we use the 
measurement of Gf^ presented in Ho et al. (2008), and the 
likelihood code described in Dossett, Ishak & Moldenhauer 
(2011). This likelihood code expands on that presented in 
Ho et al. (2008) by including the effects of modified gravi¬ 
tational held equations. 

The density held for the cross-correlation is approxi¬ 
mated by the following measurements: the 2MASS Two Mi¬ 
cron All Sky Survey, the Sloan-Digital Sky Survey Luminous 
Red Galaxy Sample, the Sloan-Digital Sky Survey Quasars, 
and the NRAO VLA Sky Survey. And the CMB temper¬ 
ature data is taken from WMAP-5^. The final Z-range we 
adopt is 6 < Z < 130: this range is taken to ensure linear 
theory is valid, specifically, this Z-range is imposed to ensure 
a wavenumber cutoff of Zs ^ 0.05 h Mpc“^. 


5 MCMC ANALYSIS 

We sample the parameter space of cosmological param¬ 
eters using Markov Chain Monte Carlo techniques with 
the CosmoMC package. The MCMC algorithm implemented 
within this code is an adaptive Metropolis-Hastings method 
which utilizes a number of techniques to ensure fast conver¬ 
gence times. The definitions and adopted priors of each pa¬ 
rameter are given in Table 2. Our results are derived using 8 
separate chains which are run until convergence is achieved. 
The convergence of the Markov chains is determined using 
the Gelman and Rubin convergence criteria, for which chains 
require R — 1 < 0.02 to be satisfied for the least-converged 
orthogonalized parameter; R being the ratio of the variance 
of the chains’ mean and the mean of the chains’ variances 
(Gelman & Rubin 1992). The posterior mean and 68% con¬ 
fidence intervals are then computed using thinned Markov 
chains. 

There is currently no consensus on the Hq value as 
measured from Cepheid data. The most up-to-date mea¬ 
surements are presented by Efstathiou (2014), Riess et al. 
(2011), and Humphreys et al. (2013): they measure Hq — 
70.6 ± 3.3, 73.8 ± 2.4, 72.0 ± 3 km/s/Mpc, respectively. Note 
both Efstathiou (2014) and Humphreys et al. (2013) have 
used the revised geometric maser distance to NGG 4258 (as 
presented in Humphreys et al. 2013), however their mea¬ 
surements still do not agree: the disagreement can be traced 
to different outlier rejection criteria being applied. For this 
analysis we adopt two approaches, because of this tension. 
When the expansion history is described by ACDM we do 
not include any Ho prior as the model-dependent constraints 
from the CMB are sufficient. When we do include deviations 
from ACDM in the expansion history we add an Hq prior 
using the measurement by Efstathiou (2014). 


^ Note, the NVSS radio survey is the best tracer of large-scale 
structure at a high-redshift: this survey provides the most signifi¬ 
cant detection of a cross-correlation. Furthermore, the ISW effect 
is only dominant at low-Z and hence is limited by cosmic variance. 
For both reasons, the measurement of has not been signif¬ 

icantly improved from Ho et al. (2008), hence justifying our use 
of this data. 


(37) 
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Table 2. Cosmological parameters used in our analysis. For each we give the symbol, uniform prior range, value taken in the ACDM 
cosmology, and summary definition. The parameters with a specified prior range are treated as free parameters in the MCMC analysis, 
while the remaining parameters are fixed at their fiducial values. The first block contains the standard parameters present in the ACDM 
model, while the second and third contain the parameters introduced to allow modifications from General Relativity. Note a prior is 
included on the derived parameter Hq. 


Parameter Prior range Baseline Definition 


. [0.005,0.1] ... Baryon density today 

Dc = Qch^ . [0.001, 0.99] ... Cold dark matter density today 

1OO0MC . [0.5,10.0] ... 100 X approximation to r»/D a 

T. [0.01, 0.8] ... Thomson scattering optical depth due to reionization 

ria . [0.9,1.1] ... Scalar spectrum index (ko = 0.05Mpc“^) 

ln{10^®Aa). [2.7, 4.0] ... Log power of the primordial curvature perturbations (/cq = 0.05 Mpc“^) 

Qk . 0 Curvature parameter today 

y~[ nil, . 0.06 The sum of neutrino masses in eV 

A^efl . 3.046 Effective number of relativistic degrees of freedom 

Parameters model I 

zt . 1.0 Transition redshift for GR modifications 

kt . 0.01 Transition wavenumber for GR modifications (Mpc“'^) 

Giigbt(fc,2) . [—10,10] 1 Modification to relativistic Poisson equation (Eq 7) 

Gmatter(fc, z). [—10,10] 1 Modification to non-relativistic Poisson equation (Eq 7) 

Parameters model II 

uiQ . [—3.0,1.0] —1 Dark energy equation of state, w{a) = wq + (1 — a)wa 

Wa . [—3, 3] 0 Redshift-dependent modification to the equation of state (see above) 

7 . [Oj 2] 0.55 Power-law index of the growth-rate parameter f{z) = 

S7a . ■ ■ • Dark energy density divided by the critical density today 

. ... Matter density today divided by the critical density 

(Tg. ... RMS matter fluctuations today in linear theory 

Ho . [20,100] ... Expansion rate today in kms“^Mpc“^ 


5.1 Parameter Fits: Model I 

Using different combinations of the measurements outlined 
in the previous sections, we performed hts to the base 
ACDM parameters (wb, Wc, 0 mc, G tis, j4s) and the modified 
gravity parameters Gmatter(fc, z) and Giight(fe, z). Recall each 
modified gravity parameter is binned in both redshift and 
scale. 

In addition to the physical parameters, a number of 
nuisance parameters are introduced to account for unknown 
astrophysical effects. For the WiggleZ multipole calculation 
for each redshift bin we include the galaxy bias and ve¬ 
locity dispersion as nuisance parameters, that is, 6 iin (2 = 
0.44), (Jv{z = 0 . 44 ),&iin (2 = 0.73), and ct„(z; = 0.73). The 
uniform priors imposed on these parameters are &iin € [0.5,3] 
and (j„ G [0,10]h“^Mpc. For the DRll-BOSS CMASS mul¬ 
tipole measurement we also include galaxy bias and velocity 
dispersion as free parameters, &iin(z = 0.53), (j„( 2 : = 0.53). 
Additionally for BOSS, we include a free parameter to ac¬ 
count for the shot noise N, this is given the prior N € 
[0, 2000]/i“^Mpc^. For the WiggleZ measurement the shot 
noise contribution has already been subtracted. For the ve¬ 
locity power spectrum measurement we include a velocity 
dispersion parameter apv(z = 0) G [0,500]km/s. 

In order to understand the sensitivity of each cosmolog¬ 
ical probe to the physical parameters, and test for residual 
systematics, we analyse different combinations of cosmolog¬ 
ical probes. The different combinations are defined and la¬ 
beled in Table 3 (henceforth we will use these definitions). 
The final results of this section are displayed in Fig. 4 and 
Fig. 5, and further information is provided in Table 4. The 


first figure shows the constraints on Gmatter(fe, 2 ) and the 
second on Giight(fe, 2 ). The black-dashed lines in both fig¬ 
ures show the predictions from General Relativity. We do not 
plot the 2D contours between Gmatter(fc, 2 ) and Giight(fc, 2 ) 
as their correlations are small, i.e., (jpcj) ~ 0.15. Here pc 
is the cross-correlation coefficient, and {) indicates the av¬ 
erage over all the possible values between Gmatter(fc, 2 ) and 
Giight(fc, 2 ). Similarly, we do not plot the inferred constraints 
on the base ACDM parameters as, with two exceptions, the 
base ACDM parameters are not highly correlated with the 
post-GR parameters, the exception being ag and with 
Gniatter- When averaging over the four Gmatter parameters 
we find (jpcj) ~ 0.77,0.39, respectively. The remainder of 
this section will involve a discussion of the content of these 
plots, in addition to some comments on potential systemat¬ 
ics effects and the derived astrophysical parameters. 


As shown in Fig. 4, we observe very little variation in 
Giight(fc, 2 ) as we add extra datasets to the base sample (the 
green contour): this is because the ISW effect on the T-T 
power spectrum is dominating the fit; additionally, galaxy 
velocities have no sensitivity to Giight(fc, 2 ), so we expect the 
benefit of including them to be minimal. The grey contours 
in Fig. 4 are derived by adding the T-g measurements to the 
base sample, and the red contours are derived by adding the 
multipole and velocity measurements to the base sample. 
And the blue contours show the main results which are de¬ 
rived using Set 4. From these measurements for Giight(fc, 2 ) 
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<l;k <0.01) Giight(^<l;^>0-01) <0.01) Giight(^ >1;^ >0.01) 


Figure 4. 68% and 95% confidence regions for the four Giight(^ 5 ^) bin parameters. Here 2 > 1 is referring to the redshift range 
2 > 2 > 1. Note all of the parameters specified in Table 2 are being varied in this analysis, however for clarity we only plot the 
constraints on Giight(^? this plot. Recall we have defined Base as High—/ + low—/ + WP + BAO + SNe. 


we infer (in terms of 68% CLs) 


G'iight(2 > 1; fc > 0.01) 
Giight(2 < 1; fc < 0.01) 
Giight(2 < 1; fc > 0.01) 
Giight(2 > 1; fc < 0.01) 


1 0 ^ 7+0053 
i.UO/_o.045 5 

1.048 ±0.048, 

1 1 CQ + OOSO 

i.iOO_0.068 5 
1.016 ±0.026, 


These measurements are compatible at the 95% CL with 
GR. 

For Gmatter(fc, 2 ) we observe a significant amount of 
variation as new measurements are added to the base sam¬ 
ple. In Fig. 5 the green, grey, red and blue contours corre¬ 
spond respectively to measurements using the dataset com¬ 
binations Set 1, Set 2, Set 3, Set 4 (sets 6 to 9 are used 
for systematics checks to be discussed in the next section). 


As derived from Set 4 (i.e. using all the datasets) the ID 
marginalised results for Gmatter (in terms of 68% CLs) are 


Gmatter(a < 1; fc > 0.01) 
Gniatter(2 < l;fc < 0.01) 
Gmatter(2 > 1; fc > 0.01) 
Gmatter(2 > 1; fc < 0.01) 


0.65 ±0.43, 


J^.^^-0.34 I 

0.53 ±0.32, 
0.87 ±0.30. 


Similarly to above, these results are consistent with GR 
at the 95% CL, while at the 68% CL level we observe a ten¬ 
sion with GR in the high-redshift and large-wavenumber bin. 
Furthermore, the constraints from Set 4 on the 2D CLs of 
the low-z high-fc and high -2 high-fe bins of Gmatter show a 
tension with the standard model at greater than 2cr. For the 
ID marginalised results this tension is significantly reduced 
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Figure 5. 68% and 95% confidence regions for the four Gmatter(fc, 2 :) bin parameters. Here 2 > 1 here is referring to the redshift range 
2 > 2 > 1. Note all of the parameters specified in Table 2 are being varied in this analysis yet for clarity we only plot the constraints on 
Gmatter: (k, z). Recall we have defined Base as High—/ + low—/ + WP + BAO + SNe. 


as the high -2 and low -2 Gmatter bins are highly correlated, 
as can be seen in Fig. 5. This degeneracy occurs as some 
probes, such as the CMB, are sensitive to integrated quanti¬ 
ties over redshift, such that higher growth at high -2 can be 
compensated for by lower growth at low- 2 . 

Introducing direct PV measurements the constraints 
shift from the green to the grey contours. The most promi¬ 
nent shift occurs in the low -2 and low-fc Gmatter bin, as ex¬ 
pected: we find a shift from Gmatter(2 < l\k < 0 . 01 ) = 
to Gmatter(2 < 1; k < 0 . 01 ) = 1.321^;^^. We 
find further improvements in the constraints for the high- 
wavenumber and low-redshift bin. Future PV surveys should 
be able to considerably improve on this situation (cf. Koda 
et al. 2014). Using the best-fit parameters from Set 4, we 
measure XmfgSv = 778 with 979 data points: the full 


6 dfGSv velocity field is smoothed onto a grid with 979 non¬ 
empty elements (cf. Johnson et al. 2014). 

Including RSD measurements results in the shift from 
the grey to red contours, for which we find a significant im¬ 
provement in the constraint on the high -2 and high-fe Gmatter 
bin. Moreover, we find that the RSD measurements have 
more influence on the high -2 bin than the the low -2 bin: 
this is an further consequence of measuring integrated quan¬ 
tities. As a systematic check we isolate the measurements 
from WiggleZ and BOSS and perform separate fits, we find 
that the two separate constraints on Gmatter are consistent. 
We can also assess how well our model fits the observations. 
By adding the multipole likelihoods we find Ax^ = 322, for 
a total of 324 measurement points. Individually, for the fit 
to the WiggleZ multipoles, with 126 data points per redshift 
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Table 4. Cosmological parameter constraints for Model I. The constraints are derived from four different groups of cosmological probes, 
we labels these groups Set 1 to 4 and define each in Table 3. For each parameter in each group we provide the 68% confidence levels. To 
keep the table a reasonable size we only consider the parameters most relevant to our analysis. 



SET 1 

SET 2 

SET 3 

SET 4 

Parameter 

68% limits 

68% limits 

68% limits 

68% limits 

Gmatter(^ < 1; fc > 0.01) 

0-96lo.L 

0.481°;®® 

0.66 ±0.47 

0.65 ± 0.43 

Gmatter(2 < 1; fc < 0.01) 

n 81 “i“0- 

-| qr) + 0.42 
-^•'^^-0.29 

-| nr) + 0.41 
-‘-•oz_0.30 

-1 Qrv + 0.39 
-*-•^^-0.34 

Ginatter(^ > 1; fc > 0.01) 

-1 OO+0.71 
-^•^'^-0.28 

-1 1 q+0.81 

-‘-•-L^_0.33 

0.54 ± 0.35 

0.53 ±0.32 

Gmatter(^ > 1; fc <C 0.01) 

Q qc:+0.42 

0.88 ±0.37 

0.82 ±0.32 

0.87 ± 0.30 

Giight(2 > 1; fc > 0.01) 

l.UO/_o 046 

1.UOD_q Q45 

1 n7’o+0.063 

043 

1 nir'7+0.053 
l.UO/'_Q Q45 

Giight(2 < 1; fc < 0.01) 

1.048 ±0.048 

1.044 ±0.050 

1.048 ±0.048 

1.048 ±0.048 

Giight(2 < 1; fc > 0.01) 

-1 1 Q+0.10 

-*-•-*^^-0.078 

-| 1 -| q + 0.098 

084 

1 14 + 010 
l-l^_0.077 

1 -icq+O.OSO 

i.lOO_Q Qgg 

Giight(2 > l;fc < 0.01) 

1.015 ±0.026 

1.016 ±0.027 

1.016 ±0.026 

1.016 ±0.026 

. 

0.02228 ± 0.00025 

0.02227 ± 0.00025 

0.02226 ± 0.00025 

0.02230 ± 0.00025 

. 

0.1172 ± 0.0013 

0.1172 ±0.0013 

0.1168 ±0.0013 

0.1163 ±0.0013 

T. 

0.087 ± 0.013 

0.0891°;°}® 

0.089 ±0.013 

0.086 ±0.012 

ln(10i°Aa) .... 

3.076 ± 0.026 

3.081 ± 0.025 

3.080 ±0.025 

3.073 ±0.025 

Ha. 

0.7011 ± 0.0077 

0.7009 ± 0.0078 

0.7031 ± 0.0074 

0.7060 ± 0.0074 

. 

0-8 . 

0.2989 ±0.0077 
0.85ltO;J^i 

0.2991 ± 0.0078 

n ■700+0.095 
u. i oo_Q Qg 3 

0.2969 ± 0.0074 

n 71 '7+0-018 
^•^l'_ 0.022 

0.2940 ±0.0074 

0 711 +9.017 
^ -‘--‘^- 0.020 

Ho . 

68.49 ±0.61 

68.46 ± 0.63 

68.61 ± 0.59 

68.83 ±0.61 


Table 3. The dataset combinations we use for fits to Model I, in 
addition to the labels we adopt to refer to them. The correspond¬ 
ing datasets should be clear from the information given in Table 
1. We define Base as the combination High—1 -|- low—Z -|- WP + 
BAO + SNe. Below CMASS refers to the monopole and quadrupole 
multipole measurements from the BOSS-CMASS sample. And 
WiggleZ refers to the monopole, quadrupole and hexadecapole 
measurements from WiggleZ (as presented above). 


Label 

Description 

SET 1 

Base 

SET 2 

Base + Direct PV 

SET 3 

Base + CMASS (fcmax = O.lOZi Mpc“^) 
WiggleZ (fc = 0.15fc Mpc“^) -|- Direct PV 

SET 4 

Base + CMASS (fcmax = O.lOZi Mpc”^) + 
WiggleZ (fc = 0.15fc Mpc“^) -|- Direct PV 
+ ISW-Density -|- CMB Lensing 

SET 5 

Base + ISW-Density 

SET 6 

Base + CMASS (fcmax = O.lOh Mpc“^) 

SET 7 

Base + CMASS (fcmax = 0.15Zi Mpc“^) 

SET 8 

Base + WiggleZ (fcmax = 0.15/i Mpc“^) 

SET 9 

Base + WiggleZ (fcmax = 0.19h Mpc“^) 


bin we measure XwiggieZ = 129.88 for the low -2 region, and 
XwiggieZ = 121.6 for the high -2 region. Finally, for BOSS, 
given we are fitting to fcmax = O.lO/i Mpc“^ , there are 72 
measurement points and we find Xcmass = 72.6. 

5.1.1 Astrophysical Parameters and Systematic Cheeks 

When calculating the power spectrum multipole predictions, 
we assumed a linear bias factor and linear perturbation the¬ 
ory. The validity of both assumptions may be questioned. 


We examine, albeit crudely, the importance of these as¬ 
sumptions by determining the sensitivity of the parameter 
fits to the small-scale cut-off fcmax. For our model fits us¬ 
ing the CMASS and WiggleZ multipole likelihood calcula¬ 
tions we ran new Markov chains using different cut-off val¬ 
ues fcSax'^®® = 0.10,0.15/1 Mpc-^and = 0.15,0.19/i 

Mpc“^. The results showed no statistically significant shift 
when the fitting range was changed. 

The astrophysical parameters for the multipole and di¬ 
rect PV fits only vary slightly when using different dataset 
combinations, hence we choose to only present results from 
Set 4 (given in terms of 68 % CLs). For the fit to the WiggleZ 
multipole we find CTv(« = 0.73) = 2.30l} g /i“^Mpc, (Jv (2 = 
0.44) = 4.468l);g/i"^Mpc, 61(2 = 0.44) = 1.089 ± 0.042, 
and 61(2 = 0.73) = 1.207 ± 0.059. For the fit to CMASS 
we find crv (2 = 0.57) = 61(2 = 0.57) = 

2.055 ±0.084, and iV(Shot Noise) = 705 ± 200/i“®Mpc^ Fi¬ 
nally, from the fit to the velocity power spectrum we deter¬ 
mine the 95% upper limit a-py{z = 0) < 334.6km/s. With 
different fcmax values adopted, one should not necessarily 
compare our results for the shot noise and velocity disper¬ 
sion with previous analysis; however, we find our bias mea¬ 
surements to be consistent with previous analysis. 

5.2 Previous Measurements: Summary and 
Comparisons 

Below we briefly summarise recent work in this field, with a 
focus on results that adopt a similar parameterisation. 

• Daniel & Linder (2010) presented constraints on {Q,V}, 
our {Giight, Gmatter}, iu bius of time and wavenumber. To 
constrain these parameters they used the following probes: 
WMAP7, supernova Union2, CFHTLS weak lensing data, 
temperature-galaxy cross correlation, and the galaxy power 
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spectrum. They identify the CFHTLS survey as responsible 
for a 2a tension with GR in the high-fc and low-2 bin for V. 
This featnre is not observed when using the COSMOS data 
or in subsequent analysis of the final CFHTLenS catalogue. 
Note RSD information was not included, therefore the final 
constraints on V are of order ~ 1. 

• Simpson et al. (2013) measure the parameters {S,p} (i.e., 
{Giight, Gmatter}) usiug tomographic weak lensing measure¬ 
ments from CFHTLenS and RSD measurements of fag from 
GdFGS and WiggleZ, in addition to WMAP7 (including low- 

1) and geometric information (see also Dossett et al. (2015) 
and Zhao et al. (2012)). Their measurements are consistent 
with GR: they find fi — 1.05 ±0.25 and E = 1.00 ±0.14. For 
this fit they assumed E, p are scale-independent and adopt a 
specific functional form for their temporal evolution: this ef¬ 
fectively confines deviations to very low-redshifts. Measure¬ 
ments of the t-g cross correlations, CMB-lensing, and the 
growth rate measurement from CMASS were not included 
in these fits. 

• Planck Collaboration et al. (2015) have recently provided 
the state-of-the-art measurements of post-GR parameters, 
placing constraints on a extensive range of specific and phe¬ 
nomenological models. For the phenomenological model they 
adopt the parameters {fj,,ri}, as implemented in MGCAMB. 
Motivated by f{R) models, a specific functional form for 
the redshift and scale dependence of these parameters is 
assumed. As appropriate to their aim, they ensure their an¬ 
gular cuts to the tomographic shear-shear measurements 
from CFHTLenS isolate the linear signal (see their Fig 

2) . This approach is not adopted throughout, however. 
Their adopted fag{z = 0.57) measurement (by Samushia 
et al. (2014)) was derived by fitting to the monopole and 
quadruple of the correlation function on scales larger than 
25h,“^Mpc. As highlighted by the authors (see their Fig. 7) 
non-linear terms are significant on these length scales, the 
result is a dependence on non-linear physics. 

In relation to the most up-to-date measurements, our 
results can be distinguished in two main ways: Firstly, the in¬ 
clusion of the velocity power spectrum measurements, which 
improve low-fc constraints; secondly, the methodology we use 
to analyse RSD measurements, and the range of RSD mea¬ 
surements analysed. We argue that the methodology of di¬ 
rectly analysing the power spectrum multipoles allows con¬ 
straints to be derived that are more widely applicable to 
non-standard cosmological models. This is because it al¬ 
lows one to restrict the analysis to scales within the linear 
regime, where the phenomenological model we use describe 
physical models (see Sect. 1). Moreover, the multipoles con¬ 
tain scale-dependent information, which is necessary if scale- 
dependent terms are introduced. 

5.3 Parameter Fits: Model II 

We now explore fits to a new parameter space that is more 
rigid regarding the allowed deviation to the growth his¬ 
tory. Two scenarios will be considered when fitting for these 
parameters, firstly, an expansion history fixed to ACDM; 
and secondly, an expansion history than can deviate from 
ACDM via. a time-dependent equation of state. We define 
the two parameter spaces as pi = {7, Wb, Wc, ^mc, t, ris, As}, 
and p2 = {7, w,iUa,cUb,a;c,6MC,r, ria. As}. We choose not 


to include the influence that deviations in the expansion 
history have on the expected growth (that is, the relation 
7 = /(70, Wo, Wa)) as the corrections are currently small. 

Note that by changing the growth rate we modify erg, 
this effect is included by altering the growth history well into 
the matter dominated regime. The modified growth factor 
is calculated as 

D{aef!) = exp f— f daHm(a)°'®®/aV (39) 

V / 

now we scale the fiducial prediction crf“*(2high) to find the 
modified amplitude crg(2eff): 

^8"(^eff) = . (40) 

The first set of results, which assume a ACDM expan¬ 
sion history are shown in Fig. 6. This plot shows the 68% 
and 95% 2D likelihood contours for the parameter combina¬ 
tions {H,n, 7} and {r, 7}. The expected value of 7 from GR is 
given by the grey-dashed line. In addition to the growth rate 
and AP constraints, these measurements are inferred using 
high-Z ± WP ± low-2 BAO (which we label in this section 
as base). For fits in this section we do not use the low-Z CMB 
T-T data or CMB lensing, since we have not included the 
dependence of these signals on 7. For the final constraint we 
measure 7 = 0.665 ± 0.0669, which is consistent with GR at 
the 95% C.L. 

The results for p2 are presented in Fig. 7, where 
we plot the 2D likelihood contours (68% and 95%), 
and the marginalised ID probability distributions for 
7, uiQ, Wa, Hm, T. Again, the black-dashed lines indicate the 
values expected from the standard model; namely, 7 = 0.55, 
wq = —1, and Wa = 0. The degraded constraint on 7 is a 
direct result of the degeneracy between the expansion and 
growth histories: this is the reason we consider both a fixed 
and non-fixed expansion history. 

We use four different dataset combinations to constrain 
these parameters, they are defined as follows: fit 1 is the 
base sample, fit 2 is base ± WP, fit 3 is base ± WP ± 
SNLS, and fit 4 is base ± WP ± JLA. We define base here 
as the combination High-Z ± Ho ± RSD/AP ± low-2 BAO. 
We use two SN samples in order to understand how sensitive 
the growth index is to our choice of adopted dataset. 

We will first discuss the main results, which are found 
using fit 3 and 4 (the red and blue contours in Fig. 7) and 
then consider how the constraints are influenced by the dif¬ 
ferent probes. Using fit 3 we infer (in terms of 68% CL) the 
marginalised constraints 

wo = -0.98l);/J, (41) 

Wa = -0.42to f7, (42) 

which are consistent with the standard model. In terms of 
deviation to the growth history we measure (in terms of 68% 
CL) 

7 = 0.761°:°*? • (43) 

This result is at tension with GR at a level greater than 2a. 
Changing our SN sample to the JLA sample we find this 
tension is slightly reduced. Using fit 4 we now measure 

7 = 0.73l?'°* , (44) 

which is just consistent at the 2a level, and for the expan- 
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Figure 6. 2D marginalized posterior distributions for and {t, 7 }, assuming a ACDM expansion history. The contours are the 

68 % and 95% CL. Below we refer to base as the dataset combination high-Z + WP + low-z BAO. The green contours are found using 
base + CMASS, the grey contours are found using base + WiggleZ, the red contours are founds using base + 6 dFGS, and the blue 
contours show the combined fit to all the growth rate measurements plus the base measurements. Moreover, we include the AP and BAO 
information with the growth rate constraints, without double counting BAO measurements. 


sion history we find wo = — 0 . 891 q'J ^2 — “0-63lo'.45- 

Note, without including any SN data, using fit 2 (the grey 
contour in Fig. 7) we measure 7 = 0.69to'ii, which is con¬ 
sistent at the 95% C.L. This may suggest there exists a mild 
tension between the growth rate and the SN measurements. 
Finally, we note our measurements of the growth index are 
relatively insensitive to the polarization data, as can be ob¬ 
served in Figure 7 by comparing the green (no WP) and 
grey (including WP) contours. 

Comparing the best-fit values for the expansion his¬ 
tory using only BAO measurements with the BAO -I- SN fit 
(which is driven by SN) is interesting as it provides a test of 
the significance of non-linear structure on SNe distance mea¬ 
surements (Clarkson et al. 2012). With fit 2, which only uses 
the low-redshift BAO measurements to constrain the expan¬ 
sion history, we infer (in terms of 68 % C.L) wo = —0.68i!1q' 26 
and Wa = These measurements are consistent 

at the 95% CL with the standard model and the constraints 
from the SN -|- BAO fit; moreover, they highlight the current 
necessity of type la SN in placing tight constraints on the 
redshift evolution of the equation of state. By introducing 
the Lyman-a BAO measurements into this fit we measure 
Wq = — 0.581q'22 and Wa = — l- 55 to'g 9 , which indicates a 
tension with the standard model predictions at a level > 2a, 
in agreement with the results by Font-Ribera et al. (2014). 
Further checks for systematics will be required to confirm 
this result given its significance and the complexity of the 
measurement. 


5.4 Comparison with Previous Results 

Below we summarise a subsample of previous measurements 
of the parameters {wo,'U’a, 7 }. 


• Beutler et al. (2014) measure 7 = 0.772lQ'Jgy using the 
power spectrum multipoles from the DRll CMASS sample 
and Planck: this fit includes the AP effect, but does not 
allow for deviation in the expansion history. This value is 
consistent with the measurement by Sanchez et al. (2013) of 
7 = 0.64±0.26 found using the clustering wedges of CMASS 
combined with BAO and SNe measurements. 

• Rapetti et al. (2013) perform fits to {uio, 7 } and 7 . For a 
fixed expansion history, using WMAP combined with galaxy 
cluster data from ROSAT and Chandra, they measure 7 = 
0 . 4151 q'J^ 26 ' When adding further data from RSD measure¬ 
ments (WiggleZ and 6 dFGS) they find 7 = 0.570lg;o63. 

• Beutler et al. (2012) measure 7 = 0.547 ± 0.088 using 
WMAP7 and the two-point correlation function measured 
from 6 dFGS. For this fit the expansion history is fixed, as 
the AP effect is not relevant. Note, there is a small difference 
between our measurement of 7 from 6 dFGS and this result. 
This change is driven by the preference for a higher flm in 
Planck compared to WMAP. 

For this analysis we extend the range of RSD measure¬ 
ments used to constrain 7 relative to Sanchez et al. (2013); 
Beutler et al. (2014) and Rapetti et al. (2013). Moreover, 
relative to Rapetti et al. (2013) we also use the updated 
Planck measurements as opposed to WMAP. The final accu¬ 
racy of our measurement of the growth index improves upon 
Sanchez et al. (2013) and Beutler et al. (2014), given the ad¬ 
ditional measurements we analyze. Note, our constraint on 
the growth index disagrees with Rapetti et al. (2013) as we 
use different datasets, and the two measurements have simi¬ 
lar accuracy as we choose to focus only on growth rate mea¬ 
surements from RSD: we do not include additional probes 
sensitive to the growth rate. This position is motivated by 
recent suggestions that there exists some tension between 
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Figure 7. 68 % and 95% confidence regions for the most relevant parameters describing model II. The base sample of datasets, as refereed 
to above, represents the combination High-1 -|- Hq + RSD/AP + low-z BAO. 


the predictions from a Planck cosmology and RSD measure¬ 
ments (e.g. Macaulay, Wehus & Eriksen 2013). 


6 CONCLUSIONS AND DISCUSSION 

In search of departures from the standard cosmological 
model and clues towards possible extensions, we have mea¬ 
sured time- and scale-dependent deviations to the gravita¬ 
tional field equations of General Relativity. We model these 
deviations using the time and scale-dependent parameters 
{Gmatter, Giight}. These parameters are defined using 2 bins 
in time and 2 bins in scale. Gmatter modifies the gravita¬ 
tional interaction for non-relativistic particles, and hence 
alters structure formation, while Gught acts equivalently 
for relativistic particles, thus affecting how light propagates 
through the universe. 

To measure the eight parameters describing this model. 


plus the six describing the standard model, we utilize a range 
of cosmological probes including BAOs, Type la SNe, the 
CMB, CMB lensing, and the cross-correlation of the CMB 
with large-scale structure probes. In addition, we include 
measurements of the power spectrum multipoles from the 
WiggleZ and CMASS galaxy redshift samples, and the veloc¬ 
ity power spectrum from 6dFGSv. Our motivation for adopt¬ 
ing a phenomenological model is to provide a set of results 
that can self-consistently be used to test the widest possible 
range of models. To this end, we have focused on only ana¬ 
lyzing measurements on scales within the linear regime. We 
summarise our main results as follows: 

• We perform a new measurement of the power spectrum 
multipoles of the WiggleZ survey, featuring a new calcu¬ 
lation of the window function convolution effects and an 
improved determination of the covariance from N-body sim¬ 
ulations. 























Searching for Modified Gravity 19 


• Modeling deviation from General Relativity in terms of the 
growth of large-scale structure, we find the following results, 
given in terms of 68% CLs: Gmatter(2 < 1; fc > 0.01) = 0.65± 
0.43, Gn,atter(2 < Rfc < 0.01) = 1.22l°®®, Gnaatter(2 > 
l;fc > 0.01) = 0.53 ± 0.32, Gn,auer(2 > l;k < 0.01) = 
0.87 ± 0.30. These constraints are consistent with GR (i.e., 
Gmatter = 1) at the 95% confidence level. We observe a small 
tension (> Icr) for the high-wavenumber and high-redshift 
bin. 

• Modeling deviation from General Relativity in terms 

of light propagation, we derive the following constraints, 
given in terms of 68% GLs: Giight(2 > l;fc > 0.01) = 
1.0571° );“, Giight(2 < l-,k < 0.01) = 1.048 ± 

0.048, Giight(2 < l;fc > 0.01) = 1.1531° “«°, Giight(2 > 
1; fc < 0.01) = 1.016 ± 0.026. These constraints are con¬ 
sistent with General Relativity at the 95% confidence level: 
the significant improvement in constraining power, relative 
to Gmatter, is due to the the sensitivity of the ISW effect and 
GMB lensing to deviations in Gught. 

• Adopting an alternative model, we introduce deviation 
in the expansion and growth histories simultaneously by 
varying the growth index and two parameters describing a 
redshift-dependent equation of state. For this ht we utilize, 
among other probes, recent growth rate constraints from 
RSDs, as measured from the WiggleZ, CMASS, and 6dF sur¬ 
veys. Our final result assuming a AGDM expansion history 
(in terms of 68% GL) is 7 = 0.665 ± 0.0669, while allowing 
the expansion history to deviate from AGDM we measure 
7 = 0.691° °°. Both these results are consistent with the 
standard model; however, introducing SN measurements to 
this ht (either SNLS or JLA) we hnd a ~ 2 (t tension with 
AGDM. 

Probes of the velocity held of galaxies have an indis¬ 
pensable role to play in addressing questions of the nature of 
dark energy as they are uniquely sensitivity to only temporal 
perturbations. The observational datasets we have analyzed 
are consistent with a vacuum energy interpretation of dark 
energy; however, due to the magnitude of current uncertain¬ 
ties any hnal conclusions drawn from these, and other cur¬ 
rent, observations would be premature. In future analysis 
tomographic weak lensing and galaxy-galaxy lensing mea¬ 
surements will be included to improve our constraints; fur¬ 
thermore, we will begin assessing the viability of specihc 
models using the inferred parameter constraints. 
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